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In this talk I review my recent work with Brian Greene and Thomas Rooslil. First 
I discuss the effect of a negative cross-coupUng on the inflaton decay in two scalar 
field theories. Our main finding is a new effect, the negative coupling instability, 
which leads to explosive particle production and very fast inflaton decay. Then I 
discuss the consequences of this instability for grandunified baryogenesis models. 
The novel aspect of this review is an intuitive explanation of the negative coupling 
instability, using field trajectories in the configuration space. 



1 Introduction 

The problem of the inflaton decay and the subsequent Universe reheat has 
received a lot of attention in recent years. At this conference an extensive 
discussion has been devoted to this problem (see talks by D. Boyanovsky and 
S. Yu. Khlebnikov) § 

A very brief history of the problem of the inflaton decay is as follows. In 
1990 Traschen and BrandenbergerH discovered that the inflaton, when it starts 
oscillating after inflation, decays exponentially fast via parametric resonance. 
The inflaton decays in a few dozens oscillations, which is_typically many orders 
of magnitude faster than the perturbative decay rate U. An example is the 
inflaton decay via a quartic coupling g(fp'x^ /2 to a second scalar field x- The 
tree-level decay rate is simply 

nH^xx) = ^ — , (1) 

where $ is the inflaton amplitude, and lo^ its frequency. Since typically g<i> <S 
cij0, the decay time Tdocay ^ is many orders of magnitude greater than the 
expansion time ^/H ^ l/ujfj,. 



"Daniel Boyanovsky presented a new inflationary scenario in which the classical limit 
gives an accurate description of the problem. Sergei Khlebnikov argued that generically, in 
theories with weak couplings, the classical limit is reached when "occupation numbers" of 
the field modes become large. 
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2 Parametric resonance in a toy model: the Mathieu equation 

Consider the simple two field case with the potential: 

mx) = ^m|02 + fxV. (2) 

In a static universe the linearized mode equations for x can be written as the 
Mathieu equation 

^ + [A{k)-2qcos{2z)]xk = Q, (3) 

where z = uj^t, and — k'^ + is the frequency squared of Xk- The 

instability chart of the Mathieu equation is depicted in figure ^. The = 
curves divide the chart into stable and unstable regions. For g > 0, A > 2\q\, 
and ^max < 0.3, as can be seen in figure |^. One can define ^ as follows: the 
amplitude of the corresponding mode grows in one oscillation of the inflaton (jj 
by a factor exp 27r/x. Since the mode occupation numbers grow as the amplitude 
squared, this growth can be interpreted as particle production, and is termed 
paramptric resonance. On the other hand for 5 < 0, A > — 2|g|, and one can 
showla that fj. < A\q\^^^/Tr, which can be 3> 1 when \q\ ^ 1. This phenomenon 
is what we mean by the negative coupling instability. The question is can this 
instability occur in a natural physical setting? 

3 Parametric resonance in a realistic model 

To understand this question, consider the general renormalizable two scalar 
field potential with a, Z2 x Z2 global symmetry 



y (0, X) = lml<f>' + \mlx' + ^^' + ^x' + f 0^X^ (5) 



in an expandiiig universe. The results I discuss in this section are obtained 
numerically in Id. Assume a natural hierarchy of couplings (1) A^^ > g > 
(when g > 0), and (2) A^ > \g\ > A^, r = X^X^/g"^ > 1 (when g < 0). In the 
former case the inflaton rolls down in inflation and then oscillates along the 
X = direction, leading to the standard growth of the resonant modes (in the 
X direction). As a consequence, the inflaton amplitude $ decays exponentially 
fast into the x fluctuations with the decay time of order 

Tdccay ^ : /i < 0.3 (6) 

2^uj^ no 
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where nmax '^1/5 and no ~ 1/2 are the maximum and initial occupation 
numbers, respectively, and lu^ is the frequency of (f> at the end of inflation. 

When g < 0, the valley of the potential (^) is not any more along the 
X = 0, but along 

Xl = {^^ -^+#g<0 . (7) 

1 otherwise 

If xo ~ Xo initially, then it remains so during inflation (as long as \q\ ^ 1/10). 
More generally, the valley (j^) is the attractor for the zero mode. However, 
during inflation the superhorizon modes grow (to generate seeds for structure 
formation, just like in ordinary inflationary theories). On flgure ^ we show 
the valleys of the potential both for = and > 0. In the former case 
the valley trajectories are defined by Xo(i)/0o(O = tan9 = ±(— f;/Aj^)^/^, and 
one can easily show that along the valley directions the initial stages of the 
evolution look just like the positive g case with q^s = 2\q\ and ficS £ 0.3. 
Indeed, we observe this behavior in our numerical simulations. However, at 
certain point the backreaction effects from created particles become important: 
the X excitations become massive with m^^g = m"^ + 3A^(x^) + g{(t>^), and 
the picture changes dramatically. The valley equation obeys (0) so that for 
= \g\4>{tY the valley has a characteristic curved shape as seen in figure |[ 
If at this time (j){t) is growing, and x has a small oscillating amplitude whose 
phase matches onto the decaying solution, then {0(i), x(0} does not follow the 
valley in (^ but instead it climbs the ridge. Consequently, the infrared x rnodes 
" see" an inverted harmonic oscillator and their amplitude grows exponentially 
fast. This is what we mean by the negative coupling instability. In this case the 
mode amplitudes and field variances grow typically much faster than in the case 
of parametric resonance. Indeed, in our numerical simulations we observe the 
instability exponent < \q\^^'^ for \q\ ^ 1. This instability occurs generically 
when <I>^ > cff — kT^^^^i and hence even in the massless x case due to 
the backreaction effects. Note the stochastic nature of the instability. In order 
to predict how each of the modes behaves, one must be able to compute the 
phase of xo at \g\(jP'(t) = very accurately. This sensitivity is illustrated 

in figure ^ in which we show the dramatic change in the variance growth, 
when the x rnass changes by a tiny amount. For example, when |g| — 350, 
Am^/m^ ~ ±0.01 leads to a change in the variance growth by a factor ~ 10^. 
The instability peaks at = 5.5 x 10~^^A/| for which ^max — 20 2± 

Note that a similar instability occurs when < 0, g > 0. In this case 
an investigation of the configuration space shown in figure ^ indicates that the 
instability occurs quite generically. Indeed, if at the end of inflation (/)q ^ 
—m^/g and x is sufficiently close to the valley value at xo = 0, once (/)o(i)^ < 
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—rn^^lg the infrared x modes become unstable, leading to rapid mode growth 
and particle production. A detailed study of this effect is under investigation. 

Without going into details, we now summarize the consequences of the 
negative coupling instability: (1) the inflaton decays much faster, typically 
within a few oscillations (with (/i) ~ 0(1)); (2) as a consequence of the following 
relation A > 2\q\ + \q\^^^, which specifies when the resonance shuts down, the 
field variances grow larger than in the positive g case by a factor ~ 4151^/^, and 
hence have a much stronger symmetry restoring force; (3) massive x particles 
are produced much easier, which may be relevant for grandunified baryogenesis 
(see below) , and the unstable momenta are of a much broader range {Ak^ ~ 
\g\^^~ml). 

4 Baryogenesis and the negative coupling instability 

In this section I discuss how the negative coupling instability can be used 
to facilitate grandunified baryogenesis. I first summarize some of the results 
concerning production of heavy particles, and then, on an example of a simple 
toy model, I discuss how this heavy particle production can affect grandunified 
baryogenesis. 

For a realistic choice of couplings, which is in chaotic inflationary models 
constrained by the COBE satellite measurements to be « 3 x 10"^'^ or 

< 2 X 10^'^GeV, one can produce massive particles with ~ lO^^GeV 
(as typically required by GUT baryogenesis models), provided \g\ > 10~^ and 
A;^ > 3 X 10~* (for stability), see figure ^ which corresponds to the parameter 
space 10^^ > \g\ > 10~^, 1 > A^ > 3 x 10""', leaving plenty of opportunity 
for baryogenesis model building. This is to be contrasted with the positive g 
case for which g > 10~^ is required. Such a large value leads to an unpleasant 
fine tuning problem since the small value of A^ needs to be protected against 
radiative corrections. 

To obtain an estimate of the baryon asyminetry that could be produced 
during preheating we study a simple toy model H. In short, the model can be 
described as follows. Initially a certain amount of energy is transferred from 
the inflaton to a heavy GUT scalar field via nonpcrturbative mechanisms, an 
example of which is the negative coupling instability. The result is a cold 
(far from equilibrium) fluid of massive particles (with a mass M^), which we 
assume decays in a B and CP violating manner into light degrees of freedom 
that instantly thermalize. This out-of-equilibrium CP violating decay is the 
crucial step that leads to net baryon production. In order to treat the problem 
analytically we further assume that at all times the temperature T of the 
light relativistic particles is ^ M^. The main results of this investigation 
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are expressions for the final baryon-to-entropy ratios ng/s for a massless and 
massive inflaton. Here we quote the result for the massless inflaton decay 



where e an effective CP violation (the decay of each pair X — X produces 
e baryons), and p\ are the initial energy densities of the X and (f) fluid, 
respectively, ^ 10'^ — 10'^ is the number of relativistic degrees of freedom, 
Mx is the mass of X, F^, Tx are the decay rates for (f) and X, respectively, and 
Hq is the expansion rate. In addition we have assumed a "natural" hierarchy 
of time scales Hq^T^^T^. The remarkable feature of this result is that the 
final baryon-to-entropy ratio does not depend on the decay constants Tx and 
r^. Why this is so can be understood as follows. As long as the temperature 
never rises close to Mx so that the massive particles X are never re-populated 
by scattering off the radiation fluid, most of the massive particles decay out of 
equilibrium and create ~ eiV^ baryons, independent on the decay rate T^. The 
entropy production, if dominated by the late inflaton decay, is also independent 
on the decay rate F^, so that the final baryon-to-entropy ratio is independent 
on the decay rates. We believe that this feature lends more credibility to our 
result. 

According to (^) the final baryon to entropy ratio is proportional to p^. 
Since the energy density of the (non-relativistic) X particles is proportional to 
their variance [px ~ ^ixii^'^Y))^ ^^'^ spikes in figure ^ are directly imaged 
to spikes in baryon production. For a massless inflaton the final baryon to 
entropy ratio is expressed in terms of the maximum variance 

From figure ^ we see that the largest value for the variance {(5X'^Y')q ^ 
IQ-^Ml, where Mp ~ 2.4 x lO^^GeV is the reduced Planck mass, leading 
to a maximum value for the baryon-to-entropy ratio tlb/s ^ 10~'^e (we set 
5* — 100). We conclude that with our toy model it is not hard to obtain 
baryon-to-entropy ratio consistent with observation n^/s ~ 2 — 12 x 10^^^. 
Certainly more realistic models for GUT baryogenesis at preheating will have 
to be constructed to verify the viability of this scenario, but the preliminary 
results are very promising. 
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Figures 

1. The stability chart for the Mathieu equation. The dark regions are stable. 
In the unstable (light) regions we marked contours of constant instability index 
^l. The contours shown are ^ = 0, 0.1, 0.2, 0.3, 0.5, 1, 2, 3. The lines A = ±2\q\ 
are also plotted. 

2. The valleys of the potential (|^). In the massless case oscillations are along 
Xo(*)/</'o(*) = tanO = ±(— g/A^)^/^, while in the massive case the trajectories 
are given by (|^). 

3. The X field variances for three runs with slightly different x masses (m^ = 
7.2 X 10-"m2, = 10-12^ Ax = 10-^ g = -10"^) 

4. The valleys of the potential with negative to^ and positive g. In the massless 
case oscillations are along xo — 0; while in the massive case the trajectories 
are given by (|^). 
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5. Maximum peak and valley variances as a function of \g\ for m-^ = lO^^GeV. 
The inflaton is massless (m^ = 0), = 3 x 10~^^, g < and A^^ is adjusted 
to keep r = 10. 
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q=350, r=10, m^^=5.3x10 Mp', <(5x) > 
q=350, r=10, m\=5.4x10'^^Mp^, <(5%)^> 
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q=350, r=10, m^^=5.5x10"' Mp', <(8x) > 
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Figure 6. Valleys of the potential, m^<0, g>0" 
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r=10,^^=3x10-'^<(6x) >peak/Mp 
r=10, ?i^=3x10-^^ < (6%)%a||ey/Mp2 
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